Influence of electron scatterings on thermoelectric effect 

Jing Lijj Tin Cheung Au Yeung, and Chan Hin Kam 

Microelectronic Division, School of Electric and Electronic Engineering, 
^ ■ Nanyang Technological University, Singapore 639798, Singapore 

(Date textdate; Received textdate; Revised textdate; Accepted textdate; Published textdate) 



o 

(N 



c^ ■ Abstract 

In this work, we employed non-equilibrium Green's function to investigate the electron transport 
CN ! properties in the nanowire with the presence of scatterings. The scattering mechanism is modelled 

by using the concept of Biittiker probe. The effect of electron scattering is analyzed under three 
(—] ■ conditions: absence of external field; with a bias voltage; and with a finite temperature difference. 
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It is found weak and strong scatterings strength affect the electron transport in different ways. In 
the case of weak scattering strength, electron trapping increase the electron density, hereafter boost 
the conductance significantly. Although the increment in conductance would reduce the Seebeck 
coefficient slightly, the power factor still increases. In the case of strong scattering strength, 
electron diffraction causes the redistribution of electrons, accumulation of electron at the ends of 
the wire blocks current flow; hence the conductance is reduced significantly. Although the Seebeck 
^,^ coefficient increases slightly, the power factor still decreases. The power factor is enhanced by 

^ . 6% — 18%, at the optimum scattering strength. 
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3] . The efficiency of a TE device 



I. INTRODUCTION 

Thermoelectric (TE) device is a type of green energy devices and is of great possibihty 
to be used widely in near future. Currently, the limitation of TE device is its low energy 
conversion efficiency, and has attracted much attention [1 
is indicated by the dimensionless TE figure of merit (ZT), determined by the expression 
S'^aT/n using its parameters Seebeck coefficient S, electric conductivity a, thermal conduc- 
tivity K (= Ke + Kph, contributed by electrons and phonons) and the temperature T. Among 
these parameters, S, a and Hf. are related to the electron transport and Uph is related to 
the phonon transport. Therefore, the strategies of developing high ZT device are improving 
the electronic TE efficiency and reducing the lattice thermal conduct ivity[lj]. Currently, 
most of the efforts are put in to reduce the lattice thermal conductivity by investigation 
the phonon transport, which can be found in our previous works Ill4l3l|. However, more 



attentions needs to be paid on the electron transport, as the electronic TE efficiency limits 
the performance of TE device significantly [ij]. The best TE device has a Dirac delta shaped 
electron transmission function|8|. Besides the TE efficiency, the power factor (S'^a), related 
to the power rating, is another important index of TE device. 

In nano-scaled device, the concept of ballistic transport of electrons is important. How- 
ever, electron scatterings still occurs for the reason of defects, grain boundaries, phonons, 
and charges in the device. Electron scatterings affect the transport properties, and further 
influence the TE effect of the device. A previous work has studied the effect of nano-particle 
scattering on TE power factor |9|. It is found that the electron concentration is usually higher 
in the sample with nano-particles, implying the Seebeck coefficient is usually unchanged and 
conductivity is increased at the optimum of the power factor. That work is based on low 
nano-particle concentration, and that work rise our interest in the problem of the electron 
scattering strength on TE effect. In this work, we employed non-equilibrium Green's func- 
tion to investigate the electron transport properties in nanowire with the presence of electron 
scatterings. Scattering points are inserted into the nanowire, and are modelled by Biittiker 



probes[lO|. The response of the device to an external field, due to a temperature difference 



or a bias voltage, is discussed with a wide range of scattering strength. 
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FIG. 1: The schematic of the modelled system: a nanowire coupled by two contacts at both ends, a 
scattering point denoted by the square box is coupled to Biittiker probe through a virtual contact. 

II. THEORY 



A silicon nanowire with square cross-sectional shape coupled by two contacts, which 
are served as thermal and electric reservoirs, is considered as shown in Fig.l. The wire is 
composed of many slices with the thickness (a^) 0.3nm parallel to the contact interface. For 
the slice at position z, the effective potential is given by 
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The variables ^^1,(2) is the quantum number for the energy in transverse direction; lx,{y) ^^^ 
side length of the cross-section (set to 5nm in the computation); V{z) the electric potential 
energy; m* the effective mass of silicon; and h the Dirac constant. The Hamiltonian of the 
wire for the mode, denoted by the subscript (ni,n2), is constructed as follow: 
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where tz is the coupling energy between two shces determined by the expression /i^/(2m*a^) 
(about 461meV in this model). 

To reveal the effect of the electron scattering, the scattering points modelled by Biittiker 
probes are inserted into the wire. At a scattering point, the slice of the wire is coupled to 
the Biittiker probe with the coupling strength (?7), as shown in Fig.l. The coupling strength 
determines the electron scatterings strength. For the convenience in discussion, we define 



the scattering factor a, the relative scattering strength respected to the couphng strength 
between two shces. Hence U = a ■ tz, and a equal to 0.056 is equivalent to the thermal 
energy at 300K, which is about 26meV. 

All the contacts in the system, including both virtual contacts (Biittiker probes) and two 
real contacts, are labeled by index i. The self-energy due to the contact i is given by 

_pi-ki(E)-az 

UE) = U' ■ , (3) 

where ki{E) is the longitudinal wave vector determined by 



h{E) = ^ ^[E - E^,,r.M]- (4) 

The broadening function of the contact i is then obtained from its self-energy: 

r,(E) = -2Im[S,(E)]. (5) 

The total self-energy of the wire S is the sum of the contribution from each contact. The 
retarded and advanced Green's function of the wire, G^ and G^, are then determined by: 

G^{E) = [E + ie-H-E{E)]-\ (6) 

G-^(E) = [G«(i?)]t, (7) 

where e is an infinitesimal positive number. With that, the electron transmission function 
between contact i and contact j is given by 

T,,,(E) = Tr[r,G'%G'^]. (8) 

The electric current of contact i is determined by the Landauer-Biittiker formula: 

I^ = ^ldEj2T,,.{f,-f,), (9) 

where fi is the Fermi-Dirac distribution at contact i with the Fermi level Ef^i, defined by 

/. = [e^ + 1]-^ (10) 

To fulfill the charge conservation, the net current of each Biittiker probe is set to zero by 
adjusting the Fermi levels of Biittiker probes using Jacobian method, (linear temperature 
gradient is assumed in the wire). The Fermi levels in the left and right contacts are deter- 
mined by assuming all donors are ionized, and the donor concentration lO^^cm"^ is used in 



the computation. The number of electrons in each shces of the wire are obtained from the 
diagonal elements of ra, determined by the following expression: 
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The potential energy V{z) for each slice is updated by the Poisson equation using the electron 
density determined by the above formula, and the final status of the system is obtained after 
the self-consistent iteration process. 

III. DISCUSSION 

In the following, the system at room temperature in three conditions: (A) absence of 
external field, (B) with a bias voltage, and (C) with a finite temperature difference, are 
discussed. 

A. Absence of external field: electron trapping and diffraction 

In this section, the length of the nanowire is set to 18nm, and two scattering points are 
inserted at the positions oi z = 6nm and z = 12nm. With the absence of external field, i.e. 
both temperatures and voltages at the two contacts are equal, the electron charge density 
along the wire determined by Green's function is shown in Fig. 2, in which the scattering 
points are denoted by "SPl" and "SP2". In the case of ballistic transport (a = 0), elec- 
tron scattering only occurs at the contact interfaces, but not in the wire. For this reason, 
electrons are accumulated at both ends of the wire, that incur the repulsion of electron as 
the result of the Coulomb force. Therefore, a fluctuation of electron density is observed in 
the wire, as shown by the "Ballistic" curve in Fig. 2a. When the scattering strength getting 
slightly stronger {a < 0.1), the electron wave function in the wire is reconstructed due to 
the perturbation of the scattering points, but this effect is not significant. Electron scatter- 
ing occurs when electron collide with the scattering points, and electrons are accumulated 
in the wire. Fig. 2a shows the stronger scattering strength, the more electrons are trapped 
■nside the w.e. especially near the seattet.ng points in the weak scattering strength region 
(a < 0.1). This result is consistent with the previous work[9]. However, in strong scat- 
tering strength region, this effect vanished. This is because the electron wave function is 
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FIG. 2: The electron density in a nanowire with the length ISnm, two scattering points are inserted 
at z = 6nm and z = 12nm, (a) is for the weak scattering strength region and (b) is for the strong 
scattering strength region, a is the relative scattering strength. 

significantly modified by the scattering points. The probability of finding an electron near 
the scattering points become much smaller, therefore the electron density near scattering 
points reduces (Fig. 2b), the chance of electron scatterings decreases, the chance of electron 
scatterings decreases. This phenomenon is known as the electron diffraction, electrons prefer 
avoiding the scattering points than being scattered. 

B. A bias voltage 



In the following, we discussed a 9nm nanowire at room temperature with a bias voltage. 
A single scattering point is considered at the center of the wire. In the case of ballistic 
transport (a = 0), due to the external field, electron density is tilting toward one end of the 
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FIG. 3: A 9niTi nanowire with the presence of bias voltage AT/, (a): the electron density, and 
(b): the electric current in weak scattering strength region; (c): the electron density, and (d): the 
electric current in strong scattering strength region, a is the relative scattering strength. 



wire, as shown in Fig.Sa. With slightly increases in the scattering strength, more electrons 
are accumulated in the wire, especially near the scattering point. This is consistent with 
the discussion in section A. With the bias voltage, the increment in electron density is 
conducive to boost the electric current as shown by Fig.Sb. However, in strong scattering 
strength region, electrons are bypassing the scattering point by reconstructing their wave 
functions. This leads to the accumulation of electrons near the contact interfaces (Fig. 3c), 
which block current flow. So, decrease in electric current is expected (see Fig.3d). 

With different bias voltages, the electric current respect to the scattering strength are 
in the similar profile: a rise followed by a drop (Fig.4). Also, the scattering strength at 
which the maximum current occurs, increases with the bias voltage. This phenomenon can 
be understood in the following way. The electrons are trapped in the wire in the weak 
scattering strength region. With the increase in the bias voltage, these electrons can be 
accelerated and contribute to the electric current. Higher bise voltage has stronger power 
in accelerating electron, therefore the scattering strength at which the maximum current 
occurs is increasing with the bias voltage. With sufficiently large scattering strength, the 
diffraction of electron causes the accumulation of electrons at both ends of the wire, which 
blocks the electric current. 
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FIG. 4: A 9nm nanowire with various bias voltages, the electric current is plotted with respect to 
the relative scattering strength a. 

C. A finite temperature diff'erence 

To discussion the TE effect, a 9nm nanowire under a temperature difference is considered. 
The temperature difference is spht into half and applied to the contacts, i.e. T^ = T + AT/2 
and Tfi = T — AT/2. In general, TE effect is evaluated in the open circuit condition. 
Therefore, the net current is zero. For this reason, to offset the current produced by the 
temperature difference, an electric field is generated. The ratio of the voltage difference 
over the temperature difference gives the Seebeck coefficient. To determine the conductance 
of the wire in the presence of temperature difference, a testing current (a small current) 
is injected into the wire. From the change of the voltage difference, the conductance is 
obtained. 

In the weak scattering strength region (a < 0.1), Seebeck coefficient is smaller with 
stronger scattering. That is because the electric conductance increases (Fig. 5b), and a 
smaller electric field is sufficient enough to offset the current produced by temperature 
difference for case of larger conductance. Similarly, in the strong scattering strength region, 
Seebeck coefficient increases as the result of the decrease in the electric conductance, as 
shown in Fig. 5c and Fig.5d. 

In Fig.5b, the conductance increases with the temperature difference, and starts to de- 
crease after the critical temperature difference, also the critical temperature difference in- 
creases with the scattering strength. This is because, when temperature difference is small. 
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FIG. 5: A 9nm nanowire with temperature difference AT, (a): Seebeck coefficient, and (b): elec- 
tric conductance in weak scattering strength region; (c): Seebeck coefficient, and (d): electric 
conductance in strong scattering strength region, a. is the relative scattering strength. 



the increment in temperature difference intends to inject more electrons into the wire; these 
electrons are trapped easily by the scattering point, lead to the a higher electron density 
and a larger conductance. When the temperature difference is more than the critical value, 
more electrons are injected than the wire could conduct, that causes the redistribution of 
electrons, which limits current flow. Due to this effect on conductance, Seebeck coefficient 
decreases at first and start to increase after the critical temperature difference (Fig. 5a). In 
the strong scattering strength region, the fluctuations in conductance and Seebeck coefficient 
are not as much as that in weak scattering strength region. This is because the chance of 
electron scatterings is reduced due to the effect of electron diffraction, number of electrons 
trapped in the wire is smaller than the case of weak scattering strength. 

Fig. 6 shows the TE power factors with the change of scattering strength. The power factor 
has maximum value at about a = 0.1, and it drops dramatically after that. The explanation 
of the change in power factor can be ascribed to the change in the the conductance and 
Seebeck coefficient. At the optimum power factor {a = 0.1), comparing to the ballistic 
case, for AT equals to OK, 20K and 50K, the Seebeck coefficient drops by 8%, 13% and 
11% (Fig. 5a); but the conductance increases by 25%, 57% and 35% (Fig.5b), therefore the 
power factor is increased by 6%, 18% and 7%, respectively. Due to the shifting of the critical 
temperature difference mentioned in previous paragraph, the power factor has the maximum 
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FIG. 6: A 9nm nanowire with temperature difference AT, the thermoelectric power factor is plotted 
with respect to the relative scattering strength a with (a): small temperature difference, and (b): 
large temperature difference. 

enhancement when the temperature difference is about 20K. 

IV. CONCLUSION 

Weak and strong scatterings strength affect the electron transport in different manners. In 
the case of weak scattering strength, electron trapping increase the electron density, boost 
the conductance significantly. The increment in conductance would reduce the Seebeck 
coefficient slightly, but overall the power factor increases. In the case of strong scattering 
strength, electron diffraction causes the redistribution of electrons, accumulation of electrons 
at the ends of the wire block current flow, hence the conductance is suppressed significantly. 
Due to this effect, the Seebeck coefficient increases slightly, but the power factor decreases 
severely. 
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